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Uncertainty of the atmospheric neutrino fluxes 

M.Honda a 

institute for Cosmic Ray Research, University of Tokyo, 
Tanashi 3-2-1, Tokyo 188 JAPAN 

The uncertainty in the calculation of atmospheric neutrino fluxes is studied. The absolute value of atmospheric 
neutrino fluxes is sensitive to variation of the primary cosmic ray flux model and/or the interaction model. 
However, the ratios between different kind of neutrinos stay almost unchanged with these variations. It is unlikely 
that the anomalous ratio {u^ / v e ) b a / {v^i / Ve) mc reported by Kamiokande and Super Kamiokande is caused by 
the uncertainty of predicted atmospheric neutrino fluxes. 



1. Introduction 

After the discovery of an anomalous ra- 
tio of Ue/Vfi in the atmospheric neutrinos by 
Kamiokande, refined calculations of atmospheric 
neutrino fluxes have been made intensively by 
several authors @|]f§[@f§. The differences be- 
tween different authors have been compared in 
Ref. Q. The major differences are in the pri- 
mary cosmic ray flux and in the hadronic inter- 
action model. These are also the main sources 
of uncertainty in the calculation of atmospheric 
neutrino fluxes. 

There have been many measurements of the 
primary cosmic ray flux, but they do not agree 
with each other. We have to consider this dis- 
agreement as the uncertainty of the primary cos- 
mic ray flux. There also have been many experi- 
mental studies of hadronic interactions. However, 
the most important piece of information, the en- 
ergy spectrum of secondary particles in the pro- 
jectile region, is not well known enough for a ac- 
curate calculation of atmospheric neutrino fluxes. 
This is partly because, as high energy physics ex- 
periments moved to the colliders, this measure- 
ment has become more difficult to make. 

Adding to the above, the density structure of 
the atmosphere is a potential source of uncer- 
tainty. In many calculations, the US standard is 
often used as the 'standard' model. However, it 
is also known that the atmospheric density struc- 
ture has a latitude dependence and seasonal vari- 
ations. 



In this paper, we study how these uncer- 
tainties affect the atmospheric neutrino calcu- 
lation in the energy range of <,10 GeV based 
on Ref. D(HKKM). Note we use the one- 
dimensional approximation as HKKM through- 
out this paper. 



2. Variation of calculation Model 

2.1. Primary cosmic ray flux 

In Fig. [j] are shown the observations of cosmic 
ray protons at < 10 4 GeV by different groups. 
We consider 3 flux models in this study: high, 
mid, and low shown in the figure for solar min. 
The flux models for primary cosmic rays agree in 
<,5 GeV region, where the differences are small 
except for these caused by solar modulation. 

Cosmic rays above a few 100 MeV originate in 
Galactic space. When they enter the solar sphere, 
they are pushed back by the solar wind. As this 
effect is more pronounced for lower energy cosmic 
rays, the energy spectrum of low energy cosmic 
rays varies with the strength of the solar wind, or 
with the solar activity. However, this modulation 
is expected to be around 5% from the minimum to 
the maximum of solar activity at 10 GeV. Above 
this energy, the effect of solar activity on the cos- 
mic ray flux is very small. 

The variation due to solar activity in the 
<,10 GeV region agrees well with the expectation 
based on the Mid-flux model and conventional so- 
lar modulation formulae. The exception in Fig. |l| 
is the crosses which are the base data for the high 
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Figure 1. Observed cosmic ray proton 

flux. Crosses from Ref . , open upward triangles 
stand for LEAP|], open squares for MASS|], 
open downward triangles for IMAX[[lO|], open ver- 
tical diamonds for CAPRIC E p/M , and open hor- 
izontal diamonds for BESS IlljT pluses, closed 
squares, closed vertical diamond, closed upward 
triangles, and closed downward triangles are from 
Refs [22], [26], [25], [34], [35], [36], and [38] of HKKM 
respectively. 



flux model. We consider that the uncertainty in 
primary cosmic ray flux is small for <,10 GeV. 

Above 10 GeV, there are rather large differ- 
ences among the different experiments. This may 
be caused by the inherent difficulties in making 
the measurements. Particularly, the problems in 
the estimation of instrumental efficiency and lim- 
ited exposure time in balloon experiments. At 
this moment we have to assume that there are 
large uncertainties in the primary cosmic ray flux 
measurements in the >,10 GeV region. How- 
ever, it is interesting that recent experimental val- 
ues with super conductive magnetic spectrome- 
ters distribute rather around the Low flux model. 

It is noted that around 20% of nucleons are 
carried by nuclei heavier than the proton. We do 
not show the details but for the heavier nucleus 



cosmic rays we use the flux value of HKKM. 

Even using the one dimensional approximation, 
the geomagnetic field affects the calculation of at- 
mospheric neutrino fluxes through the primary 
comsic ray flux due to the geomagnetic cutoff. 
Principally, however, the geomagnetic cutoff is 
not a source of uncertainty since it can be cal- 
culated to very good accuracy from the measured 
geomagnetic field. 

2.2. Hadronic interaction 

In the calculation of the atmospheric neutrino 
fluxes, hadronic interaction Monte Carlo code is 
one of the most important components. However, 
there are uncertainties in the hadronic interaction 
model, due to the lack of suitable experimental 
data. The interaction models used by different 
calculations are slightly different to each other. 

A comparison of the secondary particle momen- 
tum spectrum is shown in Fig. || and in Fig. || 
for several Monte Carlo codes which have been 
used in the calculation of atmospheric neutrino 
fluxes. (This figure is taken from Ref. Q). 
The variable x used in the figure is defined as 
x = P2ndary/ Pinddent- It should be noted that 
the difference between the histograms (BGSQ 
and HKHM§) and the smooth curve (BN@) is 
not so large as it appears. Since each value is 
multiplied by the variable x, their differences are 
emphasized in the low x region. 

It is generally difficult to modify an established 
hadronic interaction Monte Carlo code due to the 
energy momentum conservation and to the corre- 
lation of secondary particles. However, in the case 
of the atmospheric neutrinos, the correlations are 
not important. We use an 'inclusive' hadronic in- 
teraction Monte Carlo code for the study of the 
variation of hadronic interactions. 

The calculation of atmospheric neutrino fluxes 
with the inclusive interaction code is similar to 
the analytic calculation. In the case of the ana- 
lytic calculation, however, it is difficult to treat 
the competition process between decay and in- 
teraction for hadronic particles. This is the most 
crucial point in the analytic calculation. There- 
fore, we make the Monte Carlo study with the 
inclusive hadronic interaction code in this paper. 

Starting with the hadronic interaction model of 
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Figure 2. Secondary 7r + energy spectrum calcu- 
lated using the interaction model used by various 
different authors (taken from Ref. ||). 



FRITIOF 1.6 and Jetset 6.3, we consider 'mod- 
erate' variations from that. The variation of sec- 
ondary particle energy spectrum is created as fol- 
lows: denoting a 'starting' energy spectrum of 
secondary particles as 




= f P ti{x) , 



where we take x = '&2nd/^inc and E is the kinetic 
energy. Moderate variation to this secondary par- 
ticle energy spectrum can be made with a param- 
eter a as 

^ = A(a)f ptl (x 1+a ) , 

where A(a) is a factor introduced for energy con- 
servation as: 

= fxf pt i(x)dx 
[a> f xf ptl {x 1+a )dx ' 

In other words, we request that the energy sum 
for each kind of particle is conserved. This modi- 
fication method works efficiently to vary the slope 
of the secondary particle spectrum at x « 0.3. 




Figure 3. Secondary tt~ energy spectrum calcu- 
lated using the interaction model used by various 
different authors (taken from Ref. 0). 



This modification also changes the number of 
particles for each kind of secondary particle (= 
/ xfpu(x)dx). The variation of a = ±0.1 roughly 
corresponds to ±10 % variation at x = 0.3, and 
to a ±20 % variation of the multiplicity. The 
variation of a from —0.2 to ±0.2 is considered in 
this paper. The variation limit of a is set as 2 
times of the possible uncertainty of multiplicity 
(±10 %) in the hadronic interaction. 

Examples of this modification are shown in 
Fig. | and Fig. | for p+air -> tt ± at E p = 32 GeV 
with the secondary particle spectrum of BN at 
~ 20 GeV from Ref. Q for comparison. One may 
notice that the secondary particle spectrum of BN 
is outside of our variation limit, even considering 
the difference of incident particle energy. To in- 
clude BN secondary energy spectrum, probably a 
different starting point is needed. 

In addition, variations of inelastic cross- 
sections from —20 % to ±20 % and fc/7r-ratio 
from —40 % to ±40 % are also considered. Note 
that we have assumed the same a for all kinds 
of secondary particles, although f p u{x) is differ- 
ent for each particle, and Oi ne i varies at the same 
ratio for all hadronic interacting particles. The 
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Figure 4. Secondary tt + energy spectrum from 
FRITIOF 1.6 and Jetset 6.3 at 30 GeV with its 
variations for a = ±0.2. For comparison, that of 
BN at » 20 GeV is also shown. 



Figure 5. Secondary tt energy spectrum from 
FRITIOF 1.6 and Jetset 6.3 at 30 GeV, with it's 
variations for a = ±0.2. For comparison, that of 
BN at ps 20 GeV is also shown. 



variation range considered here is larger than the 
uncertainty normally considered ]l3|. 

2.3. Density structure of the atmosphere 

For the density structure of the atmosphere, 
the US-standard atmosphere model is often used. 
However, since we are interested in the effect of 
the variation of atmospheric structure, we intro- 
duce a simple model with a single scale height. 
The air density is expressed as a function of height 
as: 

Pcoi , h 
P = -r— exp(- — ) . 

For the 'standard', we take p co i = 1.231 kg/m 3 
and ho — 8.4 km, such that it agrees with the 
US-standard in global features. 

We study the effect of variations from —10 % to 
±10 % both in the scale height and in the column 
density. However, when we study the variation of 
interaction model only, we use the US-standard 
atmosphere model. 



3. Variation of neutrino fluxes 

We have summarized the variation of atmo- 
spheric neutrino fluxes for the model variations 
of primary cosmic ray flux (Fig. ^|), hadronic in- 
teraction(Fig. and atmospheric density struc- 
ture(Fig. §). 

From these figures, one can see that the varia- 
tion of the primary cosmic ray flux and the sec- 
ondary particle spectrum in the hadronic inter- 
actions have a large effect on the atmospheric 
neutrino fluxes. However, the variation of other 
hadronic interaction parameters and the density 
structure of atmosphere have far smaller effect on 
the atmospheric neutrino fluxes. 

Shown in Fig. || is the variation region of the 
ratios among different kinds of neutrinos, with all 
the variations considered here. It is seen that the 
variation of the ratios is very small. All of the 
(p e + v^)/{pn ± Vfj) variations are inside of the 
boundary of ±5 % at 1 GeV. However, the small 
statistics make the ratio variation larger in the 
>3 GeV region due to Monte Carlo method. 
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Figure 6. The variation of atmospheric neutrino 
fluxes due to the change of primary cosmic ray 
model. The ratio to that with Mid primary flux 
model (HKKM) is shown. For the flux models, 



see the text in section 2.1 



Figure 7. The variation of atmospheric neutrino 
fluxes due to the variation of interaction model. 
The ratio to the HKKM neutrino fluxes is shown. 
For the parameters used in the figure, see the text 



in section 2.2 



We have also depicted the (y e + u e )/(v^ + v^) 
ratio by HKKM shifted by a factor 1.54 = 1/0.65. 
Note that SuperKamiokande have observed the 
ratio {v^/v e ) obs /(v^/v e ) M c ~ 0.65 both for sub- 
GeV and multi-GeV regions. 

4. Summary and comments 

We have studied the effect of variation in pri- 
mary cosmic ray flux, hadronic interaction model, 
and density structure of atmosphere on the atmo- 
spheric neutrino fluxes. The variation of the pri- 
mary cosmic ray flux is directly related to the ab- 
solute value of atmospheric neutrino fluxes. The 
secondary particle spectrum in the hadronic in- 
teractions also strongly affects the absolute value 
of the fluxes. Other variations considered here do 
not have a large effect on the atmospheric neu- 
trino fluxes. 

Variation of primary cosmic ray flux and/or in- 
teraction model does not cause a large change in 
the ratio of different kinds of neutrino. We have 
noted that the secondary particle spectrum of BN 



is outside of our variation limit. This may be the 
case for the interaction model used by different 
authors. However, the ratio (p e + ^ e )/(^ + v n) 
does shows a good agreement among different au- 
thors (See the comparison in HKKM). Probably 
a variation from a different starting point would 
not give a very different answer. Thus, it is dif- 
ficult to explain the ratio {v^/ve)obs/{v^/ve)MC 
observed in Kamiokande and SuperKamiokande 
by the uncertainty in the calculation of atmo- 
spheric neutrino fluxes. 

However, when applying the calculated atmo- 
spheric neutrino fluxes to the neutrino oscilla- 
tion study, the absolute values of the fluxes be- 
come important. The determination of cosmic ray 
flux to a high accuracy and the detailed study of 
hadronic interaction at x — ^2nd/^inc ~ 0.3 are 
crucial. 

The uncertainty for the arrival direction was 
not discussed so far in this paper. We shortly 
comment on the effect of the muon bending by 
geomagnetic field and transverse momentum in 




Figure 8. The variations of atmospheric neutrino 
fluxes due to the variation of atmospheric model. 
The ratio to the 'standard' model is shown. For 
the parameters used in the figure, see the text in 
section O 



Figure 9. Variation of ratios among different 
kinds of neutrino with all the variation consid- 
ered in sections 2.1-2.3 in all the combinations. 
The dashed line shows the ratio (D e + v e )/ (p^+Vp) 
of HKKM shifted by a factor of 1.54 = (1/0.65). 



the hadronic interaction. For the muon in the 
geomagnetic field, we can calculate the average 
bending angle before decay as: 



in,, c ■ eB 



(3xl0 8 (m/sec)) 



B(T) 



eB 



2.2xlCT 6 (sec) 



106 x 10 6 (eV) 

= 0.19 x 10 4 S(T) . 

Therefore, the bending angle for the typical ge- 
omagnetic field (~ 0.5 x 10~ 4 T) is around 0.1 
radian or 5 degree. This is probably not an ur- 
gent problem at this moment. 

We note that the typical value of the trans- 
verse momentum in a hadronic interaction is ~ 
0.3 GeV. Therefore, the transverse momentum is 
unimportant for > a few GeV. 
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